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O ' Abstract 

o ■ 

■ We find an unexpected iterative structure within the two-loop five-gluon amplitude in A/" = 4 
X5 - 

D . supersymmetric Yang-Mills theory. Specifically, we show that a subset of diagrams contributing 



to the full amplitude, including a two-loop pentagon-box integral with nontrivial dependence on 
five kinematical variables, satisfies an iterative relation in terms of one-loop scalar box diagrams. 



^ . The implications of this result for the possible iterative structure of the full two-loop amplitude are 
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I. INTRODUCTION 



Scattering amplitudes in gauge theory exhibit simphcity that is not manifest from Feyn- 
man diagrams Q]. This simphcity is even more striking in the case of maximally supersym- 
metric Af = A Yang-MiUs theory (MSYM). 

One of the fascinating properties of MSYM is the possible presence of iterative struc- 
tures relating amplitudes at different orders in perturbation theory. In Anastasiou, Bern, 
Dixon and Kosower (ABDK) suggested that the planar two-loop n-gluon maximally helicity 
violating (MHV) amplitude obeys the iteration 
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where M^^\e) = A^^)(e)/A^°^ is the ratio of the L-loop amplitude (evaluated in Z) = 4 — 2e 
dimensions) to the corresponding tree- level amplitude, and 

/(^)(6)^-C(2)-C(3)6-C(4)e^ (2) 

Thisproposal was inspired by the fact that the collinear and infrared singu- 

lar 0, 0, 0] pieces of M^^^ (e) are known to satisfy an iterative relation of this form. The 
statement that ^ actually holds for the full two-loop amplitude M^^^(e) is the content of 
the ABDK conjecture, which at present has been explicitly checked only for the n = 4 gluon 
amplitude P]- A similar iterative relation at three loops was recently proven in Q], also for 
n = 4. 

The four-gluon amplitude is very special in the sense that it is a nontrivial^ function of 
a single dimensionless variable x = t/s. The n = 5 case is much more complicated due to 
the presence of several independent kinematical invariants. This is also the reason why one 
would expect a richer, less rigid structure compared to the n = 4 case. 

Although the two-loop five-gluon amplitude is not known, a conjecture for it has been 
given in [lo|. The conjecture for M^^\e) has the interesting property that contains two 
classes of terms, 

Mt\e) = vi'\e) + wP{e), (3) 

(2^ (2) 

where (e) is parity even and (e) is parity odd. This is in contrast to the four- 
gluon amplitude M[^\e), which is wholly parity even. In this paper we evaluate V^'^\e) 

^ There is a 'trivial' overall factor of (st)~^ in each term in equation 
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explicitly through 0{e~^) and numerically^ at 0(e''). As one step in this calculation, we 
present an explicit (through 0{e^^)) formula for a two- loop pentagon-box integral. To our 
knowledge this is the first appearance in the literature of a two-loop integral depending on 
five kinematical invariants. In performing this calculation we have benefited from a program 
recently developed by Czakon [11] which greatly facilitates the manipulation and numerical 
evaluation of Mellin-Barnes integrals. 

In order to check the relation for n = 5, it is necessary to know the one-loop amplitude 
M^^\e) through 0{e^). Through 0(6°), it can be written as a linear combination of one- 
loop scalar box integrals, but starting at 0{e) extra terms a ppe ar which can be expressed in 



terms of pentagon integrals in D = 6 — 2e dimensions |12l . Iisl Q| • To our knowledge, the 
contribution of these extra pieces to Mg^^ (e) has not been explicitly computed. We therefore 
decompose 

M^'\e) = V,^'\e) + wPie), W^ie) = 0{e), (4) 

where V^^\e) consists of the one-loop scalar box terms and W^^\e) contains the (currently 
unknown) extra pieces from the pentagon integrals. 

In this paper we prove the remarkable fact that the V^^^ (e) pieces alone satisfy the ABDK 
relation (P), i.e. we prove that 
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(2), 



This paper is organized as follows: In section II we define the various integrals which are 
studied in this paper and review the proposal |l^ for the two-loop give-gluon amplitude. 
In section III we provide some details of the proof of our main result (0), postponing most 
technical details to the appendix. In section IV we use double-double unitarity cuts to show 
that M^^\e) must contain the term W^'^\e), and we demonstrate the intriguing fact that 
this term has very mild IR behaviour, w!-'^\e) = 0{e^^). In section V we conclude with a 
discussion of our results and their possible implications for the iterative structure of the full 
two-loop amplitude. 



^ More specifically, it is the relation lO that we check numerically at C(e°)- After this is done, it can of course 
be turned around as in section III.B and used to construct an explicit formula for V^'^\e) through ©(e"). 
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FIG. 1: The one-loop one-mass scalar box integral. See appendix A for details. 



II. FIVE-GLUON AMPLITUDES AT ONE AND TWO LOOPS 



In this section we present formulas for the planar five-gluon amplitudes that enter into the 
main result Recall that for MHV amplitudes it is very useful to normalize loop amplitudes 
by the corresponding tree-level amplitudes. Let us denote by Mg'^'' (e) the normalized L-loop 
amplitude Af\e)/Af. We use the notation 



-{ki + kij^if. 



(6) 



The one-loop five-gluon amplitude is known to be 
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M«(e) 



\ E [W«(e)]+0(e 



cyclic 



(7) 



where J^^-* (e) is the one-loop scalar box integral shown in Figure 1 and the sum is taken over 
the five cyclic permutations of the external momenta. The missing 0{e) terms in (|7j), which 
involve D = Q — 2t pentagon integrals , are not explicitly known. We therefore define 

V^'\e)^-\ E [W«(e)] (8) 

cyclic 

with the understanding that V5*'^^(e) and Mip{e) differ starting at 0{e). This difference is 
defined to be W^^\e) = M^^^(e) - vi^\e). 

At two loops, it has been conjectured [l^ that'^ 



cyclic 



(9) 



Note that we have translated this expression into the more modern normalization conventions of 
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where the various integrals are defined in Figure 2. This formula for M^^\e) can clearly be 
written as a sum of two kinds of terms, 
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cyclic 

4 E h/'="(< 



cyclic 



(10) 

fin 



which are respectively parity even and parity odd. The results of this paper provide very 
strong evidence in support of the (e) part of this conjecture, and in section 3 we show 

(2) (2) 

that the (e) term must be present in Mg (e) as well. 



III. PROOF OF THE MAIN RESULT 

In this section we present our proof that V^^\e) and V^'^\e), defined in (jH)) and (fTUj) 
respectively, satisfy the ABDK relation (jSJ 

Vr\e) = I {Vi'\e)f + f'\e)vi'\2e) -^- + 0{e). (12) 

Our method is similar to that used in to study the four-gluon amplitude M^^\e) in 
that we are able to check (0) without the need to fully evaluate any of the loop integrals in 
terms of harmonic polylogarithm functions^. This is done by deriving the necessary identities 
between various Mellin-Barnes integrals 'under the integral sign'. 

We begin in subsection III. A by verifying (0) explicitly through 0{e^^). That this works is 
already highly nontrivial for two reasons. First of all is the fact that the ansatz © for Mg^^ (e) 
has not been completely proven, although our results obviously provide strong evidence in 
its favor. Second is the fact that, even if one assumes that the M^^\e) satisfy (0), there is no 
obvious reason why this should imply that the V^^\e) satisfy The difference between 
and © first shows up at 0{e^^) because (as we show below) W^'^\e) starts contributing to 
the left-hand side of © at this order, and W^^\e) starts contributing to the right-hand side 
of (0) also at this order. 

It seems a miracle that (jSJ holds at 0{e~^) even though so many things might have spoiled 
it. In subsection III.B we show that this miracle persists through 0{e^). It seems clear that 

^ Although it is not necessary for our proof, we nevertheless present in appendix A several explicit formulas 
in terms of polylogs, in particular for the integrals /(^^"^ and l'^''' which are new to the literature. 
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A;4 ki + k2 k^ 




FIG. 2: Cast of characters: here we define the five two-loop integrals which appear in this paper. All 
figures refer to the corresponding diagrams in a scalar field theory, i.e., a collection of propagators 
integrated over the loop momenta. The integrals l(^)'^ and are multiplied by factors which 
depend on one of the loop momenta, as indicated explicitly above. These factors are meant to 
be included in the numerators of the corresponding scalar integrals, transforming them into tensor 
integrals. See appendix A for further details. 

the 0(e°) terms in (0) could, in principle, also be checked without fully evaluating the 
integrations explicitly, but it would require a larger number of more complicated identities. 
For the purpose of this paper, which is to point out an unexpected iterative structure within 
the full amplitude, we are content to perform the final step at 0{e^) numerically, via a robust 
procedure described below. 
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A. The ABDK Relation for V^^\e) Through 0{e~^) 

Explicit formulas for the integrals and /'^^^"(e) have appeared in the literature jl^ . 

Q]. We have evaluated the new integral /'•^^'^(e) explicitly through 0{e^^). The results for 
these integrals are summarized in appendix A, where they are expressed in terms of several 
basic Mellin-Barnes-type integrals which we call A(a), Aj{a) and F{a,b). The A functions 
are defined by the integrals 

dz 



A{a) = f-^a'r\-z)r{z)r{l + z), 
J zm 



2TTi 

dz 



(13) 
(14) 
(15) 
(16) 



A,{a) = f a^T\-z)T{z)T{l + zmz) +^), 
A,{a) = f^a^T\-z)T{z)Til + zm-z)+^), 
A,ia) = [^a^T\-znz)Til + z)-, 

J ZTll Z 

where il){x) = ^ Inr(x) and in each case z is integrated from —ioo to +ioo along a contour 
which intersects the real axis between —1 and 0. The function F{a,b) is given by a more 
complicated double integral defined in appendix A. Here we only need to use the fact that 
it is symmetric in a and b. Explicit formulas for all of these functions in terms of harmonic 
polylogarithms are given in appendix A. 

If we plug the expressions given in appendix A for the various integrals into (0), we find 
that the 0{e~^) through 0{e-^) pieces cancel easily, leaving just 
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+ C(e° 



, Si-2 , 



(17) 



Interestingly, all of the F(a, 6) terms automatically drop out of this expression due to the 
symmetry F{a,b) = F{b,a). The remaining terms on the last two lines of (fT7j) can be seen 
to vanish with the help of the identities 

A{a) + A{l/a) 



2 ^ ' 

Ai{a) - A2{a) + ^As{a) = -^ln(a)A(a 



TT 

y 



A,{l/a)-A2{l/a) + -A^{a) 



1 TT^ 1 

- ln(a)A(a) - — ln(a) - - \n^{a), 



(18) 
(19) 

(20) 



which can be derived by applying various tricks directly to the definitions ()13|) - ()16|) . This 
illustrates our point that it is possible to verify (fTTj) without the need to fully evaluate loop 
integrals in terms of harmonic polylogarithms. 

An interesting fact, addressed in section V below, is that the integral I^"^^^ does not seem 
to appear in V^'^\e) (in the chosen basis of integrals). Suppose we make a more general 
ansatz 

- ^ E [faS3sll^'^''{e) + hslsj^^^\e) + hs,s^s,I^'^%e) + hs,s,sj^'^\e)\ (21) 

cyclic 

and ask whether there exist numbers fa,-''yfd such that (j2II) satisfies At order C(e~^) 
and 0{e~^) the only constraint is that the /'s should satisfy 

fa + fb + Ifc + 3/d = 5. (22) 

However, at C(e^^) the functional dependence of the various integrals on the Si is much more 
complicated, and it is easy to see that the unique solution enabling (|7T|) to satisfy (0) even 
through 0(e~^) is 

fa = fb = fd= 1, fc = 0. (23) 



B. The ABDK Relation at 0{e^) and The Finite Remainder 

Having demonstrated that V^'^\e) satisfies the ABDK relation through 0{e^^), let us now 
turn our attention to the 0{e^) piece. In principle one could continue as in the previous 
subsection by identifying the basic integrals which emerge from the Mellin-Barnes represen- 
tation and then working out the necessary identities relating them to each other. However, 
as mentioned above, it suffices for our purpose here to perform this final step numerically. 



A very efficient program 
been developed by Czakon 



br numerical evaluation of Mellin-Barnes integrals has recently 



nl (see also H). ch term in ^ 

separately, for randomly generated values of the kinematical variables Sj, and always found 
that although each term in (jSj) is generically 0{1), the left-hand and the sum of the terms 
on the right-hand side conspired to agree with each other to within the precision of the 
numerical integrations (~ 10"^). This highly nontrivial conspiracy leaves us with no doubt 
that (0) is correct. 
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A robust way to distill the information contained in at 0(6") is to study the so-called 
finite remainder of ^^^•'(e). The finite remainder of an arbitrary L-loop n-gluon amplitude 
was defined in For the case at hand, we construct the quantities 

5 



= -iE(^^)"^ (24) 

/(^)(e) ^ -l(/«(e))%/(^)(e)/«(2e), (25) 
in terms of which the finite remainders for V^^\e) are defined by 



FPie) = \/W(6)-/«(6), (26) 
(e) = Vi'\e) - [P\e) + n'HeWi'^e)] . (27) 

Of course 0] defined this quantity with the full amplitudes M^^\e) instead of just the pieces 
V^^\e) we are using here. Nevertheless, the analysis of that paper proves that if V^'^\e) 
satisfies (jSI), then (e) must be given by 

Ft\e) = I iFt\e)Y + f^He)Ft\2e) - ^ + 0{e). (28) 



Of course it is simple to calculate F^^\e) from the definitions (jHI) and (j26j) . 

H\) = - J E >" f ^) 1" (^) + ^ + (29) 

4^ \Si+iJ \Si+2j 8 

Plugging this result into (|28p leads to the following sharp prediction: If (0) holds at 0(e°), 
then the quantity defined in (j2Z|) must be equal to 

where we have defined 

X = Eln^f^V F = Eln^f^). (31) 

Having checked the ABDK relation (0) as described above, the formula ()30|) for the finite 
remainder must also be true as a consequence. Nevertheless we found it useful to also 
check (jHUj) directly via a more robust numerical procedure as follows. Supposing we didn't 
know ()30|) but wanted to check it numerically, we began with the ansatz that F^'^\o) should 
be a quadratic polynomial in In^ of kinematical invariants. Now the right-hand side of (fTTj) is 
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manifestly invariant under cyclic permutations Si Sj+i, so (e) must have this symmetry 
as well. The objects X and Y defined in ()31|) are the only independent quantities with the 
right structure satisfying this symmetry. We therefore made the ansatz that F^'^\(]) should 
be a quadratic polynomial in X and Y with six coefficients to be determined, 

fP (0) = aiTT^ + aavr^X + agvr V + a^X^ + a^XY + ae^^ (32) 

The factors of vr^ were chosen conveniently so that all of the were expected to be rational 
numbers. Indeed, we found that the rational numbers oi = 89/1152, 02 = —03 = 11/96 
and 04 = as = —a^jl = 1/32 were repeatedly obtained to a precision of 10~^ for various 
randomly generated values of the kinematical invariants, so we are quite confident that the 
finite remainder F^'^\o) is indeed given by (jHUj) . The robustness of this numerical check 
comes from the fact that we use numerical data to fix just six rational numbers. 



IV. CONSISTENCY CHECKS ON THE W^ie) PIECE OF (e) 



m 



n this section we study the structure of the proposed two-loop five-gluon amplitude given 
lot . First we show that V^'^\e) cannot possibly be equal to the full amplitude M^^\e) 
by itself. This is done by computing a double-double cut in a two-particle channel to see 
that the wP{e) piece is needed. We then show that W^'^\e) has very mild IR behaviour, 
divergingonly as 0{e^^). This fact makes the presence of W^'^\e) consistent with the known 



IR lei 0,3 behaviour of M^"^(e). 



(2), 



A. Double-Double Cut 



Unitarity cuts provide very powerful constraints on scattering amplitudes. At one loop, 
amplitudes in supersymmetric gauge theories are completely determined by their cuts |lS,l3|- 
This fact is especially powerful if the basis of master integrals (after reduction has been carried 
out) is known. In particular, at one loop in A/" = 4 SYM, the problem of computing any 



20|. 



amplitude is reduced to that of computing tree amplitudes by using quadruple cuts 

At two loops, the basis of integrals is not known. In the case of n = 4 only double-box 
scalar integrals are needed to write the full amplitude j2]- This suggests that for n = 5 the 
basis should include the double boxes J^^)'', and J*^^-*^ shown in Figure 2. In addition 
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one might expect pentagon-boxes, i.e, a pentagon joined to a box along one propagator. 
These are impossible for n = 4 but for n = 5 the integral /^^•''^ shown in Figure 2 is natural. 
One should also allow for the same classes of integrals with various tensor structures in the 
numerator. 

The proposal of 10] for M^^\e) is a linear combination of the integrals in Figure 2 with 
tensor structures taking the general form 



(33) 



Our goal in this subsection is to show that the pieces in V^ie), as defined in are not 
enough to be consistent with unitarity and that the proposal of (ij] is just right to give the 
correct cuts. 

Consider the double two-particle unitarity cut in the Si channel. Recall that unitarity cuts 
compute discontinuities across branch cuts of the amplitude. In this case, the double-double 
cut computes the discontinuity of the discontinuity of the amplitude across the branch cut 
in the Si channel. This is given by (we use the notation Pjj... = ki + kj + ■ ■ ■) 

Cll = J ^ 6^^\q')6^^\{q + P,,,r)A'-{-{q + P,,,)-XA\5\q-) 

X / ^<5W(£^)5(+)((£-Pi2)')A*^^^((-g)+,r,(-£ + Pi2)-,(g + P345)+) 

xA*=--((-£)+,l-,2-,(£-Pi2)+). (34) 



The integral over i is easy to recognize as the cut in the Si channel of the one-loop 
amplitude A^^\{-q) + ,1' ,2~ , {q + Ps^^y). Therefore^ 



^12 



12 



5W(g2)5(+)((g + P3,5)2)^t-^^_^) + ^ ^-^ (g + P 



(27r)^ 

xA*^-(-(g + P345)-,3+,4+,5+,g-) 



345 J 



X 



(4vr; 



2-e 



(35) 



{k2 



^ The strange prefactor in the last hne converts from Feynman diagrams, which are normahzed according 
to J d,-^ p / {2tt)-^ , to the normahzation (jA2|) for the scalar box diagram. 
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This expression can be further simphfied by following the steps in section 5 of where 
unitarity cuts of one-loop MHV amplitudes are studied. The first step is to use that = 
and (g + ^345)^ = to write 

A'^^'^ii-qy, 1-, 2-, (g + Pu,)-')A'^''i-iq + P345)", 3+, 4+, 5+, g") = 



V^^(l-,2-,3+,4+,5 



tr_(#,(^ - If,,)^,) ii-mS - ^12)^2 



{q-h)\q-hf {q + hfiq-hf 
li-mM - Ifi2)h) tr-(#5(^ - ^12)^3 



(36) 



(g - A;i)2(g + P^^f (g + hf{q + P,^f 

where tr_(«) = |tr((l — 75)»). 

Now we expand all terms inside the bracket in such a way that the g-dependence is only 
in propagators or in tr(75...) terms. This immediately gives the result 

n 2 



^12 



1 

12-2 



(47r; 



2-e 



r(2)a 



H+2 



HsiS2 - tr(75#i^ j3))/S?+3 + {s,s, - tr(75#3^4^5))/(2)'^ 1 . (37) 



We use an abbreviated but hopefully transparent notation here: On the first line, /j^j_|_2 
means the double-box integral J^^^" shown in Figure 2, but with ki replaced by ki+2 (so that, 
for example, the massive leg has + k^ instead of ki + k2). Moreover, all three tr(75 • ■ •) 
terms are meant to appear in the numerator of the corresponding integrand. Finally, the 
three diagrams 

j(2)a^ j(2)6 j[2)d 

appearing on the right-hand side of (j37jl should of course 
be double-double cut in the obvious way. 

Note that except for the terms with tr(75...), (jHTj) would be consistent^ with vl^'^\e). 
Luckily, the first two tr(75...) terms give no contribution at all. To see this note that the 
whole integral depends on only three independent momenta and by Lorentz invariance the 
trace must give zero. For example, the first term in (jHTf) only depends on ^5, fci, and k2- 

Now it is clear why the last term does not vanish. The reason is that the pentagon-box 
depends on four independent momenta. This shows that the (e) in M5 (e) has to be 
there in order to be consistent with unitarity. 



^ Recall from equation (1) of [3| that M,^\e) is normalized with a factor of [2e '^^ /{AttY ^] relative to 
Feynman diagrams. Taking this into account turns the prefactor in (|37|l into — i, in agreement with (|10|l . 
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(2') 

B. Infrared Behaviour of (e) 

Here we discuss the very special infrared behaviour that the parity odd term (e) 
possesses. The structure of IR singularities as poles in e was first studied at two loops by 
Catani who showed that the coefficients of the and terms in a two-loop 

amplitude M^^)(e) are given by"^ 

(6) = \ (mW(6))' + /(^)(e)M«(26) + 0{e-^). (38) 

On the other hand, in this paper we have explicitly shown that (0) 

Vi'\e) = \ (M«(e))' + /(2)(e)M«(26) + 0{e-') (39) 

(we have replaced V5*"^''(e) by m''^^ [t) on the right-hand side since the difference w'^^ [e) 
only starts contributing to this equation at 0(e^^)). The only way to avoid a contradiction 
between ^ and dSl is if the remaining part ivf ^(e) = Mf\t) - V^^\t) is 0{t-^). 
By Lorentz invariance it must be true that 

Wf\e)=ix{n444^^:)wf{e) (40) 

where (e) is a function of only the Sj kinematical invariants. To write (e) in the 

form pn|) is straightforward using Feynman parameters where the presence of the original 

only leaves a Feynman parameter in the numerator. It is easy to construct a Mellin-Barnes 

representation for vj^^^{e) (see appendix A. 2). For example, one can take the integral ()A21|) 

with an additional factor of 

r(-2 - 'ie-v - Zi)r(-3 - 2e - 2:1 - 2:3 - Z4 - 2:5 - z^) ^^^^ 
r(-l - 3e - - Zi)r(-4 - 2e - 2:1 - 2:3 - 2:4 - z^ - z^) 

in the integrand. By explicit computation using the program jll| we find the remarkable 

result that 

(e) = ^ + co + 0(e) (42) 



where for example 



Cl oc 



xr(i + + 2;2)r(-i - 2:1 - 2:2 - Z3)r(i + + ^3) 

xr(-l - 2:1 - 2:2 - 2:4)1(1 + 2:2 + 2:4)1(2 + 2:1 + 2:2 + 2:3 + 2:4). (43) 



^ This is adapting Catani's formula to the case of A/" = 4 SYM. 
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This is a rather surprising result, because one would expect a generic two- loop tensor integral 
q^q^ to diverge as e~^. This integral is special because it contains the factor q^K^ in the 
numerator, where the vector K^^^ = efj,ypa{ki + k2Yk^k^ has zero inner product with many of 
the external momenta. One has to go all the way to 0{e^^) before a term can appear with 
sufficiently complicated structure to survive being dotted with K^. 



V. CONCLUSION AND DISCUSSION 

In this paper we have proven the striking result that V^^\e) and V^'^\e) satisfy the ABDK 
relation (0), even though V^^\e) is only part of the full one- loop amplitude M^^\e), and 
(e) is only part of the full two-loop amplitude Mg . Given this result, it is easy to see that 
the full ABDK relation (0) holds if and only if the 'extra' pieces W^^\e) = M^^\e) - V^^^\e) 
are related by 

W^'\e) = vPie)W^'\e) + 0{e) = -^.{s^s.sss.s.r^' w!,'\e) + 0{e). (44) 



In [l^ it was conjectured that 

wP{e) = E [^i^^'^^(^)] • (45) 

cyclic 

We have performed a couple of consistency checks on this conjecture. First, we studied 
some unitarity cuts that explicitly show the presence of the Wf'{e) term inside Ml '{e). We 
also remarked that (jl^j) has a very mild leading IR singularity C(e^^), which is consistent 
with (gH) since wP{e) is C(e). 

Finally, it is interesting to compare the formula for M^^ (e) in (jU)) with the results obtained 
by using the recent technique of octa-cuts introduced in 2l|. The basic idea of the octa-cut 
technique is that at two loops one can cut all propagators belonging to a box, i.e. perform a 
quadruple cut j^l inside a double box or a pentagon-box and get from the Jacobian a new 
propagator. The new propagator is just right so as to produce a new effective one-loop box 
or pentagon integral. The remaining one- loop object can then be cut four times to localize 



the final loop integration momentum. In 2l|, the coefficients of all double-boxes for the 
five-gluon amplitude were computed. Perfect agreement is found with Q for the J^^^" and 
j(2)b integrals. However, there seem to be a discrepancy for the /^^•*'^ integral, which is missing 



from (jUj) but should have a nonzero coefficient according to 21 1 
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The solution to this puzzle is to realize that /(^^"^ secretly contains an 7*^^)^ integral. To 
see this, expand the numerator {q — ki)'^ = (f — ^ki ■ q; the first term cancels a propagator 
and produces precisely the I^'^^^ Inte gra" 
calculation of the I^"^^^ coefficient in 



. We then find perfect agreement with the octa-cut 

9 



This preliminary analysis shows that © satisfies several non trivial constraints. As a 
future direction, it would be very desirable to make a more detailed analysis of all unitarity 
cuts of M^^\e). It would also be interesting to compute W^^\e) in order to determine whether 
the ansatz satisfies (jH}. 

APPENDIX A: EVALUATION OF INTEGRALS 

This appendix contains many of the technical details referred to in the text. The one- and 
two- loop integrals we need to study are shown in Figures 1 and 2. Our convention is that 
each loop momentum integral comes with a normalization factor of 



with D = 4 — 2e. This is the conventional normalization for studying iteration relations and 
agrees with that used in writing (P), although it differs slightly from the normalization of 
amplitudes which should appear in the actual S"- matrix (see !2| for details). As an illustration, 
the one-loop one-mass scalar box integral of Figure 1 is defined to be 



ki + k2 ks 



We begin in subsection A.l by recording our results for the various integrals, before 
detailing the method by which they were obtained. In subsection A. 2 we present Mellin- 
Barnes representations for the new integrals I'^'^^^ and J^^)'^. Finally in subsection A. 3 we 



^ Actually, a given integral can be decomposed in two different ways but we choose as the canonical 

decomposition the one just described. 
^ Slightly more subtle is to prove that © has the correct octa-cut in all "channels" . A more detailed analysis 

which is beyond the scope of this paper shows that is the case. 
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explain how we obtained the results of A.l by expressing everything in terms of a small 
number of basic Mellin-Barnes-type integrals. 



1. Results 



In what follows we use the notation 



Li = In Si 



(A3) 
(A4) 

Alia) = - a) + i/o,o.i(l - «) " ^H,{1 - a), (A5) 

o 

A^ia) = - a) + /fi,o,i(l -a) + ifo,i,i(l - + ^o,o,i(l - «) + C(3), (A6) 



and the functions 

A{a) = - 



-ifi,i(l-a)-ifo,i(l-a)--, 

o 



TT 



A3(a) = - a) + i7i,o,i(l - a) + — - a) + 2C(3), 





(A7) 



and 



1 — a — b 



+-^0,1,1 
-^0,1 f 



h-a-b 

1 — a — 6^ 
1 — a 



1 / l-g-fe 

2^°'^^^ l(l-a)(l-6), 



^1,1,1 1 



ln(l - a) + iJo,: 



1 — a — b \ , , vr^ 
,(l-a)(l-6)j Ml-a)-yln(a) 



-Ho,i (l - ln(a) - (^1 - ln(a) + ^ In^ll - a) ln(l - 6) 
+ (a ^ 6). 



(A8) 



We use here standard conventions for harmonic polylogarithm functions (see for example the 



very useful program 



22]). It is very likely that the expression ()A8|) can be simplified using 



various harmonic polylogarithm identities. However, as we have emphasized, it is possible to 
explicitly verify ()17p without knowing the precise formula ()A8|) . It turns out that as long as 
F{a, b) is symmetric in a and b, which we have made manifest in ()A8|) . then it automatically 
drops out of (fTTj) . 

The one-loop one-mass scalar box integral J^^-* was evaluated to all orders in e in Q|. In 
our notation it is given by 



jd) 



^l-e r 2 



+ 4^) + e/S + 0(e^ 



(A9) 
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with 



r(l) 



j(l) 
-'+1 



2A + 2A {'A + ; 



S3 



S4 



-2F (X 21n r^iU, f^) - 21n f 
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C(3). (AlO) 



The two-loop integral /'■^^'^ was evaluated through 0(e°) in |l6|. In our notation it is given 



by 



with 



{S3)-'-'' is,)-'-'' 



(Sl) 



-2e 



1 I'^'l'^ I'^'l" 



+ 



(All) 



r(2)a 



(2)a 



-'-1 



rii, ilU 2^. (A + fil) - 2A, fill + A, fil) - 3A, ' 



Si Si 
Si 



S3 



S4 



+A3(^)+21n^U^ +41n^U 

,S4/ VS3/ Vsa/ Vs4/ VS4 



S3' 
Si 



S4 



S3 



25 



C(3) 



(A12) 



(A13) 



The two-loop integral J^^-*^ is clearly given simply by 



/(2)^(e)(Si,S3,S4)=/(^)"(e)(Si,S4,S3). 



(A14) 



The two-loop integral J*^^-'^ has not appeared in the literature. We evaluate this integral^° 
only through C(e~^), since that is sufficient to rule out its appearance in the ABDK relation 
(see section III. A). We find 



/(2)'=(e) 



{siS3)-'-''/Hs2)-'-'' 



[S4S5J 



-2e/3 



n /(2)c 

4,4+ ,2 ^^y^ > 



(A15) 



with 



r(2)c 



S4 



S4 



S5 



S5 



1,^2 (Si 



1 



S3\ , StT^ 



Finally we have the integral j(^)'^, which also has not appeared in the literature. We find 



j(2)d/ 



(s2S5)-^-/^(s3S4)-^-/^ 
Si 



-'-2 , -'-1 



0(e°) 



(A17) 
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MS and AV are grateful to R. Roiban for collaboration on early attempts to evaluate this integral. 
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where 



A2)d 

'-2 



-3A{-^]-2a(-^]-a(-^ 



and 



A2)d 

-'-1 



S3 



S2 



Si 



sa 



S5 



^Si 



(A18) 



-• 


^ 1 


+ 2F (|, 


S4\ 


VS3 


53/ 




Si/ 



S4 S4 

.S3/ Vsi/ VS4/ Vsi 



Sl 



-A2 

-3A, ( - 



S3 



Sl/ \S4 
■^sA , . /S4 



S4 



2^2 ( - - 2^2 

-Sl/ VS2 

''^4\ , „ < /S3 



S4 



l ^i) 

S5/ 



-S3/ \Si/ \Si/ \S2/ \S5 

+(3Li - 2L2 - 4L3 + L4 + 2L5)A f-) - i(3Li + 2L2 - 2L3 + L4 - 4^5)^ 

VS3/ 3 

-i(3Li - 4L2 + L3 - 2L4 + 2L5)A (^^) + (2L2 + L3 - L4 - 2L5)A (^^ 

+ ^(L2 + 2L3 + 2L4 - 5L,)A (^^) - ^(5L2 - 2L3 - 2L4 - L,)A (^^) 

2 o_ . 1 

27 



+ ^7r2(L2 + 2L3 + 2L4 - 5L5) + 7^(2^2 + L3 + L4 - 2L5^^ 



9 

~7j(-^2-^3 + -Z^2-^4 + -^^2-^3-^4 + L^L^) + — L2(L3 + La,)L^ + -L\L^ 

+ ^(L3L4 - - L^)L5 - + 14C(3) 



(A19) 



2. Mellin-Barnes Representations 



In order to evaluate these integrals we used the technology of Mellin-Barnes (MB) integral 
representations (see [2^1 for a thorough review). It is straightforward to develop a MB 
representation for all of the above integrals (see jl^] for a fairly general treatment). For 
example, the new integrals J*^^-*^ and J*^^-*^ admit the (highly non-unique) MB representations 



/(2)^(e) 



,2e7 



.3+2e 



U=l 



X ■ 



r(-2e) 
r(-l - e- zi- Z2) 



S5/ VSs/ VSs/ \S5 

r(i + ^2)r(-i -e-23) 



r(i-^i)r(-z2)r(-i-3e-z3)" 

xr(l + zi + Z'i)T{-l - e + zi + Z2 - Zi - Z5 - Z6)T{1 + Z4 + zq) 
xr(-2;2 + Z5 + Z6)T{-2 - 2e - z-i - Z4 - zq - zj) 

Xr(l - Zi + Z4 + -27)r(3 + 2e + Z3 + Z4 + Z;^ + Zq + Z-j) 



(A20) 
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and 



lim ( 



X 



r(-2e) 
T{-l-e 



zi)r(i + zi + z2)r(-i 



S2 



32/ \S2 
Z2 - Z3) 



Z7 



r(-2 - 3e - z/ - 2i)r(l - Z2)T{iy - z^)Y{Z + 6 + ^1 + ^2 + ^3) 

xr(i + Z3)r(2 + e + 2:1 + + Z3)r(i - ^2 + ^4 + 2^5) 

xr(-4 - 2e - - - 2:4 - - -26)r(3 + e + + ^2 + 2:3 + -25 + -Ze) 

Xr(-3 - 2e - Z/ - 2:1 - 25 - - -27)r(l + -26 + z-j) 

xr(4 + 2e + z/ + 2:1 + 2:4 + ^5 + ^6 + 2:7)- (A21) 

A few comments are in order. In each case the contours of integration for the Zi variables 
can be taken to be hnes parallel to the imaginary axis, as long as Re(2;j) are chosen such 
that the arguments of all F functions have a positive real part. In the expression ()A21|) . the 
quantity v originates as the power of the factor {q — ki)"^ in the denominator of the Feynman 
integral. Since we want this factor to end up in the numerator, we want to take u ^ —1. 
This limit must be defined by analytic continuation, since it is not possible to satisfy all of 
the constraints on the Zi contours if one takes z/ = — 1 at the start. 



3. Building Blocks 

A very efficient package for developing the e series expansion of a general Mellin-Barnes 
integral has recently been developed by Czakon Using this package, we can easily read 
off the values of the desired integrals, order by order in e, in terms of some basic building 
blocks. Specifically, we find that for all of the integrals considered in this paper, only the 
following two integral structures appear through order e~^: 

A{a) = a^T\~znzni + z), (A22) 

B{a) = I ^ a^+' F(-l - z)V{-z)V\l + z). (A23) 

The integration variable z in each of these expressions can be taken to run along a line 
parallel to the imaginary axis from — ^ — too to —\ + zoo. These integrals can be evaluated 
explicitly using standard techniques. For a > we find A{d) as written in ()A4|) . and 

S(a) = i^o,i(l-«)-y- (A24) 
19 



We have chosen to express the argument of each dilogarithm as 1 — a so that the branch cuts 
of the integrals (1X221) and (1X231) for a G (—00,0) manifestly map on to the branch cuts of 
the harmonic polylogarithms on the right-hand side. 

Although we have written explicit formulas for A{a) and B{a] for the reader's benefit, one 
of the points we want to emphasize in this paper (following [15]) is that it is more efficient 
to check ABDK-type relations 'under the integral sign.' For example, instead of evaluating 
the integrals ()A22|) and ()A23|) explicitly, it is simpler to notice that by taking z ^ —1 — z 
inside the integral we see immediately that 

A{a) = B{l/a). (A25) 

This procedure of identifying the basic Mellin-Barnes building blocks such as ()A22|) 
and ()A23|) . and then deriving various identities between them, can be carried out (with 
increasing complication, of course) at each order in e. At 0{e^^) a number of new integrals, 
including double- integrals, start to appear. Using various tricks 'under the integral sign,' 
they can all be expressed without too much difficulty in terms of the basic integrals 

Alia) = J^^a^ T'{-z)T{z)T{l + zmz) + 7), (A26) 

A2{a) = / ^ «^ T'i-z)Tiz)T{l + zm-z) + 7), (A27) 

Asia) = l^.a^ r^(-z)r(^)r(l + z)-, (A28) 
J Im z 

F{a, b) = fp:p: a'+^^^b'+^^r{~l - ^i)r(-^i)r(l + ^i)r(-l - Z2) 

xT{-Z2)T{l + Z2)T{2 + zi + Z2). (A29) 

We have tabulated the explicit evaluations of these integrals in subsection A.l, but we em- 
phasize that it is possible to verify the ABDK relation ^ through 0{e^^) without needing 
to evaluate these integrals, due to the identities (fTH)). 
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